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Abstract: We present the calculation of the electroweak corrections for squark-antisquark 
pair production at the LHC within the Minimal Supersymmetric Standard Model. Taking 
into account all possible chirality and light-flavor configurations, we evaluate the NLO 
EW corrections, which are of (D(a 2 a), as well as the subleading tree-level contributions 
of 0(a s a ) and (D(a 2 ). Numerical results are presented for several scans in the SUSY 
parameter space and relevant differential distributions are investigated. The impact of the 
electroweak corrections is nonnegligible and strongly depends on the chirality configuration 
of the produced squarks. Our analysis includes a discussion of photon-gluon initiated 
processes with a focus on the impact of the corresponding large PDF uncertainties. 
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1. Introduction 

Supersymmetry (SUSY) [1] is one of the theoretically most appealing beyond -the-Standard- 
Model (BSM) scenarios. In particular, the minimal supersymmetric extension of the Stan¬ 
dard Model (MSSM) [2-4] predicts a light scalar Higgs boson that is compatible with the 
resonance observed at the LHC experiments [5,6]. Moreover, the MSSM naturally stabi¬ 
lizes the electroweak vacuum and, assuming R-parity conservation, provides also a dark 
matter candidate. 

These features give to SUSY a central role in the BSM searches at the LHC. Within 
the MSSM the analyses of the first-run LHC data have investigated several final states 
and signatures [7,8]. All the analyses found no significant deviation from SM predictions, 
setting limits on the SUSY spectrum. Still, the second run of the LHC, at higher energies 
and higher luminosity, will probe larger regions of the MSSM parameter space that are not 
excluded at the moment by experimental searches. 

At the LHC, the production of a pair of colored SUSY particles (qq',qq'*,ggiM^) 
plays a crucial role in the context of direct searches; their cross sections are the largest 
among all the possible production mechanisms. The phenomenological relevance of this 
class of processes has triggered an extensive effort in improving the precision of their 
theoretical predictions. Not only the leading-order (LO) contributions [9-13], but also 
the complete set of next-to-leading order (NLO) QCD corrections [14-20] and most of the 
NLO electroweak (EW) corrections [21-33] are available in the literature. Approximate 
next-to-next-to-leading order (NNLO) QCD corrections to squark-antisquark production 
have been computed as well [34-37]. 

Besides fixed order calculations, the large corrections in the threshold region have been 
computed at the next-to-leading (NLL) [38-42] and next-to-next-to-leading (NNLL) [43-47] 
logarithmic accuracy. Effective field theory approaches allowed also for the resummation 
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of soft-gluon and Coulomb corrections [48-51]. Moreover, gluino and stop bound states 
have been studied [52,53] and the bound-state effects in gluino-gluino and squark-gluino 
production have also been analyzed [54-56]. 

Also the squark and gluino decay rate is known at the NLO accuracy, including both 
QCD [57-59] and EW [60,61] corrections. Phenomenological studies systematically includ¬ 
ing the NLO QCD corrections to the combined production and decays of squark-squark 
and stop-antistop pairs have been respectively performed in [18,62] and [63,64], In [65,66], 
the NLO QCD corrections to squark-squark and squark-antisquark production have been 
matched to parton showers. 

From the experimental side, the lack of a signal from direct searches has increased 
the lower bounds on the masses of the supersymmetric particles, especially for those that 
are colored. In particular, the most stringent bounds have been set on light-flavor squarks 
and gluinos. For instance, assuming a simplified MSSM scenario with only the gluino 
and degenerate light-flavor squarks production with decay into a massless neutralino, the 
gluino (light-flavor squark) mass should be heavier than 1.4 TeV (1.6 TeV). Under the same 
hypotheses, degenerate light-flavor squarks and gluinos should be heavier than 1.8 TeV [67]. 

In these experimental analyses the simulation of the signal takes into account only part 
of the aforementioned higher-order corrections. In particular signal events are obtained by 
following [68], i.e., they are produced at LO accuracy and afterwards they are globally 
reweighted via the A'-factors from the inclusive NLO+NLL QCD total cross section that 
is obtained from the codes Prospino [69] and NLL-fast [42]. The theoretical uncertainty 
affecting the signal is evaluated by varying both the parton distribution function (PDF) sets 
and the renormalization/factorization scale. It is worth to notice that this procedure does 
not account for the kinematical dependence of the higher-order contributions and for the 
significant impact of the NLO corrections to the decays of the produced particles [18,62-64], 
Moreover corrections of EW origin are completely ignored both in the evaluation of signal 
and in the estimation of the theoretical uncertainty. Since EW corrections naturally involve 
most of the MSSM particles via loop corrections, they formally depend on the specific 
scenario considered. For instance, a large mass gap between colored and uncolored particles 
may in principle enhance their effect beyond the percent level. In order to have under 
control the impact of higher-order corrections, it is necessary to explicitly check these effects 
and identify for which parameters EW corrections possibly show sizable dependencies. 

In this paper we perform this type of analyses, focussing on the hadronic associate 
production of a light-flavor squark (a squark of the first two generations) together with a 
light-flavor antisquark, 


PP 


q a $ : 


( 1 . 1 ) 


where q, q' = u,d,c,s, while a and (3 denote the chirality of the corresponding squark, 
a, f3 = L, R. We compute the EW corrections to the processes (1.1) extending and complet¬ 
ing the analysis of [26], which focused only on same-flavor, same-chirality squark-antisquark 
production, PP —>• q a q *. Moreover, we present for this process the first phenomenological 
study combining NLO QCD and NLO EW corrections. Our study mainly focuses on the 
NLO corrections to the inclusive squark-antisquark production cross section, and it can 
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be directly used in experimental analyses to estimate the theoretical uncertainty related to 
the missing EW corrections. Additionally, the impact of the EW corrections on kinematic 
distributions for undecayed squarks is considered as well. 

The outline of the paper is as follows. Section 2 describes the various partonic processes 
contributing up to 0(a 2 a) and the strategy of the calculation. The numerical impact of 
the EW and QCD corrections are discussed in Section 3, followed by our conclusions in 
Section 4. 

2. Calculation method 

In our calculation we consider and treat independently all the possible light-flavor and 
chirality configurations for squark-antisquark pairs, as denoted in eq. ( 0 ). Since we treat 
quarks of the first two generations as massless, the squark chirality eigenstates are also 
the mass eigenstates, thus each squark-antisquark pair can separately be considered as a 
physical final state. In total, the possible two chirality/mass eigenstates for each of the 
four (anti)squark flavor (it, d, c, s ) lead to eight distinguished types of (anti)squarks and 
thus to sixty-four possible q a q'p squark-antisquark pairs. 

For each one of these sixty-four final states the LO cross section is of 0(a 2 ). In our 
calculation, we take into account the complete set of NLO EW corrections, which are of 
0(a 2 a), and all the contributions originating from tree-level diagrams only. The latter 
include the LO, but also 0{a s a ) contributions, which arise from the interference of 0(a s ) 
and 0(a) amplitudes, and 0(a 2 ) contributions from O(a) squared amplitudes. With the 
term “EW corrections” we will in general indicate the sum of 0(a 2 a), 0(a s a ) and 0(a 2 ) 
contributions. 

All these contributions are independently calculated for each one of the sixty-four final 
states. Thus, the masses of the eight different (anti)squarks can be treated as completely 
non-degenerate. In our phenomenological analyses of Section [3], we will include also NLO 
QCD corrections by the help of the code Prospino [69], which is based on the calculation 
in [16]. However, in order to obtain genuine NLO QCD corrections, the usage of Prospino 
and thus the evaluation of NLO QCD corrections will be performed only for degenerate 
squark masses and only at the level of inclusive cross sections. 

In the following we describe the organization of the calculation of the EW corrections 
and we explain in detail how the different perturbative orders enter for the sixty-four final 
states. To this end, it is useful to identify classes of processes that share the same type of 
partonic initial states for the different perturbative order under consideration. Among the 
sixty-four final states we may distinguish five different classes of processes: 

• ULU^-type processes producing squarks with same flavor and chirality, 

PP^u a u* a , d a d* a , c a c* a , s a s* a . (2.1a) 

• ULUp-type processes producing squarks with same flavor and different chirality, 

PP->u a u*p, d Q d* 0 , c a c*p, s a s*p , a^/3. (2.1b) 
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Processes 


Channels 



0(a s ) 

0{a) 

0[otJ 2 a 1 / 2 ) 7 -ind. 


U L U* L - type 


u l c* * l - type 


u L d * L -type 


ulu* r and 
uls* l - type 







qa 


ip 




qa 


qp 


Table 1: Representative tree-level diagrams for the different partonic subprocesses entering 
the various squark-antisquai'k production processes. The diagrams in red enter only if 
q = q', while the diagram in blue contributes only if q and q' belong to the same SU(2) 
doublet. The quarks p q and r q are defined in eqs. (|2.4|). 


• ULC^-type processes producing two up-type or two down-type squarks with different 
flavor, 

PP^u a c*p, c a u*p, d a s*p, s a d*p. ( 2 . 1 c) 

• ULSL-type processes producing one up-type and one down-type squark of different 
families, 

PP^ru a s*p, s a u*p, d a c*p, c a d*p. ( 2 .Id) 

• ULd£-type processes producing one up-type squark and one down-type squark of 
the same family, 

PP ^ u a d*p , d a u*p, c a s*p, s a c*p. ( 2 . 1 e) 

This classification is based only on the technical aspects of the computation. In particular 

it does not consider the dependence of the partonic contributions on the chirality and the 
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flavor of the produced squarks, e.g., dnd* R production is considered to be a ULU* L -t ype 
process. Therefore members of the same class entail numerically different contributions. 
For example, all necessary matrix elements for d^d* R production can be obtained directly 
from ulu* l matrix elements by flipping the chirality and SU(2) states. In general, a global 
flip of chiralities, families and/or SU(2) states project a specific squark-antisquark pair, 
and its corresponding matrix elements, into another one of the same class. In order to 
correctly perform such a flip, the Feynman rules for all interaction vertices are written in 
the most general form allowing for arbitrary chirality and SU(2) states. Correspondingly, 
for a specific chirality or flavor configuration within a given class individual diagrams might 
vanish. The presented classification allows for a dramatic reduction of the number of matrix 
elements that have to be practically computed. 

In order to specify the structure of the EW corrections, for any partonic subprocess 
X we will denote the contribution of 0(a“a b ) to the total hadronic cross section as . 
This notation can be easily generalized to any observable and differential distribution for 
the squark-antisquark system. 

We start considering the contributions to the cross section from tree-level diagrams 
only. As already noted, these contributions can be of 0(a 2 ), 0(a s a ) and 0(a 2 ). Since 
quark-antiquark initial states can produce a squark-antisquark pair with the same flavor 
configuration via a QCD or a EW gauge boson or gaugino in the t-channel, the amplitudes 
of 0(a s ) and O(a) for qq' —> q a q'p are always present. The corresponding Feynman dia¬ 
grams are shown for each class in the first two columns of Table [I]. Consequentially, 0(a 2 ) 
and 0(a 2 ) cross sections are always nonzero for each squark-antiquark chiral and flavor 
configuration. On the contrary, 0(a s a ) cross sections can be zero, since the interference 
of 0(a s ) and O(a) amplitudes can vanish in particular classes due to the color flow or 
chirality structure. 

For this reason, it is convenient to separate the contribution and a 0 ?, 

QQ -^QotQp QQ QotQp 

from all the others in the sum over all parton channels, 


2,0 _ 2,0 , - 2,0 

PP^q a q'p qq'^qaq'p PP^qnq'p ’ 


a 


. 1,1 

P P -^Qaq'p 
0,2 


- 1,1 

a PP^q a q'p* ’ 


_ 0,2 - 0,2 

(J PP^q a q' l3 * ~ (T qq'^q a q' l3 * a PP^q a qfi ' 


( 2 . 2 ) 


In eq. (2.2) all the process-dependent contributions are included in the barred quantities 
(j“’ D _ The specific structure of <r“’ D _ for each one of the five classes is listed in 

rr^q a qp ±'f^>q a qp 

the following, 


U L U* L - type: 


a 2,0 


2,0 


PP^rq a q% u gg^qa"* 


2,0 

q"q"^q a ql 

q"+q 

a PP^q n q* x — a qq^q a q^ °gi^q a q% 


- 1,1 


a - >q a q* + 

- 1>1 


(2.3a) 


- 0,2 


a 


PP^fq a ql 


E 

q"+q 


0,2 

a q"q"^q a q* a 
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u L c* l - type: 


- 2,0 _ n 
a PP^q a q'* ~ U 

- 1 - 1 — n 

°PP^q a q!S ~ 

- 0,2 _ 0,2 

G PP^>q a q!* ~ <T pqp q r-+q a q'p 


(2.3b) 


u L d*L- type: 


- 2,0 _ n 

a PP^q a q'* ~ 

- 1,1 _ 1,1 

a PP^q a q'* ~ °qq'-+q a q'Z ) 

- 0,2 _ 0,2 

17 PP—*qaq'p ~ a r q f q ,^q a q'* 


(2.3c) 


- 2,0 _ n 
a PP^q a q'* ~ U 

u L u R and u L s* L - type: °pp^q a q>* = 0 

- 0,2 _ n 
. tT PP^q a q^ 


(2.3d) 


In the previous equations the quark p q denotes the SU(2) partner of the quark q , 


Pu = d, 


Pd = u, 


Pc = s. 


Ps = c, 


(2.4a) 


while the quark r q is defined as follows, 


r u = c, 


Td = S, 


r c = u , 


r s = d , (2.4b) 


i.e. r q and q are same-charge quarks belonging to different generations. All these contri¬ 
butions have been successfully checked at numerical level against both Madgraph [70] and 
Prospino [69]. 

The NLO EW corrections, i.e. the O(ala) contribution, constitute the original com¬ 
putation of this paper. Besides virtual corrections, this order receives contributions also 
from the real radiation of gluons, photons and (anti-)quarks. The O^a^a) corrections to 
the cross section can be written as follows, 

2,1 _ 2,1 , 2,1 2,1 
a PP^q a q'* — a qq'^q a q'p ^99'^>q<xq!p 9 + ^qq! ^q a q!£l 

, 2,1 2,1 - 2,1 (r> 

+ + a q'9^q a q'$q + G PP^q a q'$ ' 

• — 2 1 

Again, with the barred quantity dAo ~ we denote the set of initial states that are 

-f" 1 rqa C[p 

class-dependent. The term dpp^- is nonzero only for ULU* L -type processes, i.e., 


u l u* l - type: 


- 2,1 _ 2,1 , 2,1 

(J PP^rq a q'S ~ °99^9<y.q.a + <T 99^q a qi l 7 


+ E {<! 

q"^q v 


q"^q a qa °q"q"-*qaq%g 


, 2,1 , 2,1 2,1 

+ <T q"q"^q a q% 7 + °9" 9 ^qa.q* a q" + °q" 9 ^qa.q* a q" 


, ( 2 . 6 a) 


other processes: er^l, ~ = 0 . 

r r 7q a cip 


( 2 . 6 b) 
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At the same order of perturbation theory, i.e. at 0(a 2 a) in principle also NLO QCD 
corrections to the 0(a s a ) photon-induced squark-antisquark production enter. We do not 
consider these contribution here. 

The NLO EW corrections ( p.5Q depend on the full set of the MSSM parameters and 
have been evaluated by using FeynArts [71,72], FormCalc [72,73] and, for the evaluation 
of the one-loop integrals, LoopTools [72], The cancellations of the ultraviolet divergences 
involves both 0(a s ) and 0(a ) one-loop renormalization and has been extensively described 
in [26,31]. Infrared singularities (IR) are regularized by introducing a small mass for the 
photon and the gluon, while quark masses are kept as regulators for the collinear singulari¬ 
ties. In the processes considered, the IR singularities of gluonic origin are Abelian-like and 
can be safely treated by using mass regularization. Infrared and collinear singularities are 
handled by the double cut-off phase-space-slicing method [74-76], along the lines of [31]. 
The initial-state collinear singularities of gluonic (photonic) origin are factorized and ab¬ 
sorbed in the parton distribution functions by using the MS (DIS) scheme. We carefully 
checked that, on the level of inclusive cross sections and of individual distributions, all 
our numerical results do not depend on the two phase-space slicing parameters and on the 
fictitious gluon and light-flavor quark masses. 

The amplitudes for partonic processes with (anti) quark-gluon in the initial state may 
exhibit an internal gluino, neutralinos or charginos that can go on-shell. This happens 
when one of these particles is heavier than one of the (anti)squarks produced. However, 
the corresponding Breit-Wigner distributions appear only at O(a^) or 0(a s a 2 ), but not 
at 0(a 2 a). Thus, similarly to the squark-squark production case [26,31], in NLO EW 
corrections these singularities do not correspond to a physical resonance. In order to avoid 
numerical instabilities, the poles have been regularized by including the width for the 
resonant particle in the corresponding propagators. Practically, the width plays the role of 
a regulator parameter, numerical results do not depend on its value. 

Technical details of our calculation have already been extensively discussed in [26], 
where NLO EW corrections were calculated for the ULU^-type processes only. 


The NLO QCD corrections, of 0(af), depend only on the mass of the squarks and of 
the gluino and have been computed in [16]. For the phenomenological studies of section 
we evaluated them for total cross sections in the degenerate squark-mass case with 
the code Prospino [69]. For completeness we show the general structure of the NLO 
QCD contributions, using the same notation adopted in eqs. (^T|) and (2.6) for NLO EW 
corrections: 


3,0 _ 3,0 . 3,0 , 3,0 

a PP^q a q'* ~ a qq’^q a q'p °qq'~^q a q'p 9 + °qg^rq a q'^q' 

, 3,0 , - 3,0 

fJ <?' g^qaq'a q + a PP^q a q'n ' 


(2.7) 


3 0 

Also in this case, the barred quantity depends on the flavor and the chirality of 

the final states: 


u L u* L -t ype: 


- 3,0 _ 3,0 , 3,0 

a pp^q a q'* ~ ^gg^qaqa + 


+ E 

q"¥=q 


3,0 

a q"q"^q a q* a 
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( 2 . 8 a) 


, 3,0 , 3,0 , 3,0 

+ a q"q"^q a q* a g + °'q"g^q a q*q" + °q"g->q a qZq" I ’ 


other processes: Cpp ~ = 0 . 

rqctqp 


( 2 . 8 b) 


3. Numerical results 


The results presented in this paper are obtained for a hadronic center-of-mass energy of 
\/S = 13 TeV. The numerical evaluation of the hadronic cross sections has been performed 
by using the NNPDF2.3QED_mc PDF set [77], which includes the photon PDF and LO QED 
effects in the evolution of all the PDF members. Moreover, this set of positive-definite PDFs 
provides an excellent description of all hard scattering data [78] avoiding the negative cross 
sections for high-mass particle production described in [79]. 

The relevant Standard Model input parameters are taken from [80] and, consistently 
with the specific PDF set used, we set a s (mz) = 0.119 in the numerical evaluations. 
Fermions of the first two generations are considered as massless. All the MSSM masses 
and couplings are determined by eleven independent TeV-scale parameters, 

m / 1 0 , tan/3 , X t , p, M 2 , nig, A% L , A%r , M~ t , Mg 3 , M^ . (3.1) 

where, for the soft-breaking sfermion-mass parameters the relations 


-^t,L/R — Mb, L/R — ^23 ) l/r — Afg l/R ) (f ~ d, C , S ) , 

= Afp T , L = , Afj^^p = A tvf,~L = Mji, (/ = e,/i), (3.2a) 


are assumed, and the bino mass M\ and the sfermion trilinear couplings are set to 


Mi 


5 sin 2 6 \v , , 

O- Tn~ ’ 

3 cos z dw 


A t = X t - 


tan {3 ’ 


A / = 0, 


(/ = e,fj,T,u,d,c,s,b ). 

(3.2b) 


In eq. o. Mi is obtained from gaugino-mass unification at the GUT scale. The physical 
spectrum and all the Feynman rules are derived from this parameter set by using tree-level 
relations. The only exception is the physical mass of the Higgs bosons, which is computed 
with the help of FeynHiggs 2.10 [81-85]. 


In order to study both degenerate and non-degenerate scenarios, in the numerical 
discussion we will consider particular slopes of the (M^l, M^r) plane, parametrized in 
terms of a mass Mg as follows (see also Fig. |l|): 


slope S\ : 

(Afg,L,Mq jR ) 

= (1, 1) • Mg , 

slope S 2 ■ 

(Afg,L,Mq jR ) 

= (1,2)-M g -, 

slope S 3 : 

(M q -L,Mg, R ) 

= (2, 1) • Mg , 

slope S 4 : 

(Afg,L,Mg jR ) 

= (Mg, 1500 GeV) 

slope S 3 : 

(M,-,L,M g - R ) 

= (1500 GeV, Mg) 
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(3.3) 






Slopes in the (Mg : l, M^a) plane 



Figure 1: Slopes in the (M^l,-M^r) plane. 


The slope S i represents the degenerate case, whereas the other slopes correspond to 
different possible non-degenerate spectra. The motivation for considering non-degenerate 
scenarios is twofold. First, EW corrections depend on the chiralities of the produced par¬ 
ticles. Thus, a mass gap between the different-chirality squarks is expected to modify the 
impact of EW corrections on the inclusive cross section for squark-antisquark production. 
Second, squarks of different chiralities may be experimentally distinguishable by their dif¬ 
ferent decay products [62]. This is already relevant for the degenerate case 1 and plays an 
even more important role in case of large mass gaps, where the production of heavy states 
is suppressed and the experimental signature is determined by the decays of the lighter 
states. 

In order to identify the dependence of the cross sections on the chirality of the produced 
(anti-)squarks, we will consider four different classes of final states: 

• “LL” : production of a left-handed squark and a left-handed antisquark, 

• “RR” : production of a right-handed squark and a right-handed antisquark, 

• “LR+RL” : production of a squark and an antisquark with different chiralities, 

• “inch” : inclusive production of a squark-antisquark pair, 

where the inclusive production is by definition given by the sum of all the final states in 
the “LL”, “RR” and “LR+RL” classes. Depending on the class and the point of the slope 
considered, the masses of a given squark-antisquark pair can be considerably different. 

1 Experimental analyses typically assume simplified models, where all the squarks decay directly into the 
lightest neutralinos. Allowing for different decays, very different signatures would emerge from left- and 
right-handed squarks. The corresponding bounds from direct searches would be consequentially modified. 
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As a general approach, we set the factorization scale fi p and the renormalization scale jUr 
equal to the average of the mass of the produced particles, gp = Mr = ( jriq a + m q' p )/2 , 
independently for each final state. It is worth to note that each one of the classes of 
processes in eq. (2.1) can in general contribute to LL, RR and LR+RL production. The 
only exceptions are ulu* l - type processes, which do not contribute to LR+RL production, 
and ulu * r -type processes which do not contribute to LL and RR production. 

In the numerical results we will refer to any contribution either specifying explicitly 
its perturbative order (and partonic initial-state) or using the following standard notation: 


• “LO” refers to O(ag) contributions, 

• “NLO QCD” refers to the sum of O(al) and 0(o? s ) contributions, 

• “EW corr.” refers to the corrections of EW origin, i.e. to the sum of 0(a s a), 0(a 2 ) 
and O(a^a) contributions, 

• “EW corr. (no g 7 )” refers to the EW corrections without the contribution of the <77 
channel. 


Before presenting our numerical results, we want to comment on the choice and the 
motivations of presenting EW corrections with and without 0(a s a) contributions from the 
gy initial state. 

3.1 PDF uncertainties in the gluon—photon channel 

Direct production of massive particles, i.e. with masses above 1 TeV, probes PDFs at large 
Bjorken-x, where they are poorly constrained by experimental data and exhibit large un¬ 
certainties. In the case of the direct production of 1 TeV (2 TeV) squark-antisquark pairs, 
e.g., the intrinsic uncertainties on the gluon and the quark PDFs lead to an uncertainty on 
the total NLO QCD cross section of the order of 4% (17%) [79]. EW corrections originate 
in general from processes with the same initial states of those from QCD corrections. Thus, 
the uncertainty on the total cross section induced by intrinsic PDF uncertainties from EW 
corrections is expected to be a fraction of the analogue contribution from the NLO QCD 
cross section. 

An important exception in our calculation is the 0(a s a ) contribution from gluon- 
photon initial state, which appears in same-flavor same-chirality squark-antisquark pro¬ 
duction (u,LU* L -ty pe processes). This initial state contributes to the EW corrections, but 
it does not contribute to LO cross sections and NLO QCD corrections. At 0(a s a), the 
only MSSM parameters entering the cross sections in the gy channel are the masses of 
squarks and the gluino. Moreover, the photon PDF in NNPDF2.3QED is affected by large 
uncertainties [77], which have to be taken into account to correctly identify the impact 
of higher-order corrections [86,87]. For this reason, before discussing the effect of EW 
corrections, we consider the PDF uncertainties arising from 57 -induced production and 
we describe their impact on our calculation and phenomenological studies. Here, we do 
not want to perform a complete analysis on the effects of PDF uncertainties from EW 
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(a) 


(b) 
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(c) (d) 

Figure 2 : Absolute and relative contributions of the photon-gluon channel to the total 
hadronic cross section. Central values (black dash-dotted lines) and PDF uncertainties 
(blue bands) are plotted as function of a common squark mass rriq ( nig = 1500 GeV). In 
panel (a) A gy->q a q* denotes the ratio between the PDF uncertainties of the cross section 
of the 57 channel and its absolute central value. 


corrections 2 , instead we want to show how the evaluation of photon-induced production 
without PDF uncertainty can give qualitatively misleading results in the context of TeV- 
scale supersymmetric-particle production. 

Fig. n a) depicts the contribution of the 57 -channel to squark-antisquark production 
as a function of the common mass of the produced squarks, rriq (degenerate case). The 
blue band around the central value (black solid line) corresponds to the PDF 68 % c.l. 

2 This kind of study is typically not done even for SM processes. In case of discovery, it can be easily 
performed if it is necessary. 
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m A 0 

700 GeV 

tan /3 

20 

x t 

2 Mq 3 


350 GeV 

m 2 

350 GeV 

m-g 

1500 GeV 


1500 GeV 

Mq.K 

1500 GeV 

Mq 3 

500 GeV 

M h 

1000 GeV 

Ml 

500 GeV 



Table 2: Default parameters within our eleven-parameter phenomenological MSSM. 


error computed along the lines of [79]. In particular, since the photon PDF is positive by 
construction, replicas distribute in a very non-gaussian way. The 68 % c.l. error around the 
central value corresponds to the symmetric error that includes 68 of the 100 results given 
by the different PDF replica from NNPDF2.3QED. In the lower panel, we plot as function of 
rriq the value of A 97 , which is defined as the ratio of the 68 % c.l. error and the central value 
for the 57 -channel. As expected, the PDF uncertainty is very large. It ranges from 75% 
at rriq = 1 TeV to almost 100% at rriq > 4 TeV. Confronted with such a huge uncertainty 
it has to be understood that for heavy squarks the contribution from the 57 -channel is 
compatible both, with zero and the double of its central value. In a new version of NNPDF 
or in another future set of PDFs with QED evolution the size of the error could be very 
different. Thus, it is wise to keep the 57 -channel always as a separate contribution to the 
overall EW corrections. 

The importance of keeping this term separate is even more evident when its size is com¬ 
pared to that of the LO cross section. An illustrative example may be found in Fig. §(b), 
which shows the impact of the 57 channel relative to the leading order contribution, com¬ 
puted assuming nig = 1500 GeV. In the following we will focus on the productions of 
squark with masses below 2500 GeV, where the 57 channel is supposed to be smaller than 
30% of the LO cross-section. This mass range will be probed by Run-II of the LHC [ 88 ]. 

Actually, the relative contribution of the 57 channel to the overall squark-antisquark 
production cross section is diluted by the “LR+RL” processes, which do not involve the 57 
initial state. In the case of only LL or RR production, which do include the 57 initial-state 
contribution from ui,u* L -t ype processes, the impact of the photon PDF uncertainty is even 
more important, cfr. Figs, ^(c) and ||(d). With rriq = 2.5 TeV, the 57 channel can induce 
corrections from zero up to 80% of the LO results. Clearly, this means that the precise 
contribution from the 57 channel is largely unknown, but it is potentially large. Conversely, 
without taking into account PDF uncertainties, one would be tempted to claim large effects 
of order 40% from photon-induced EW corrections. In the following discussion, we will not 
show the uncertainty band of the 57 channel. However, as general rule, one has to bare in 
mind an 0 ( 100 %) uncertainty is always associated to its contribution. 

3.2 Inclusive cross sections 

For our numerical evaluation we consider the benchmark scenario defined in Table [ 2 ], 
which is a slight modification of the “light-stop” scenario defined in [89]. Starting from 
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this scenario we perform various one-dimensional scans, by varying the parameters in the 
list O- Specifically, we considered scans over M 2 , m~ gi +, tan/3 and five different slopes 
in the (M^p, M^r) plane, which have been defined in (^). Since we always consider 
l/r > 1 TeV, the soft breaking parameter M g l /r can be safely considered as equal 
to the physical masses m^L/R - 3 The regions considered in the scans are within the limits 
from the Higgs sector, i.e. they exhibit a Higgs-state close to the observed one, with cou¬ 
plings compatible with the observed rates in the Higgs search channels. The compatibility 
between the scenarios and the experimental results has been checked by using the codes 
HiggsBounds [90-93] and HiggsSignals [94], by following the procedure described in [33]. 

The results of the scans are shown in Figures ||-|]. Since the scans over 5 and tan /3 
show that the total cross section is very insensitive to their values, we did not include them 
here. Each figure collects six plots related to a particular scan. Plot (a) shows the LO 
and NLO cross section predictions, inclusive over the sixty-four squark-antisquark pairs. 
In the lower panel of plot (a) we show the size of the corresponding relative corrections 
to the LO result. NLO QCD corrections are included only for degenerate squarks. The 
plot (b) shows the contributions of the individual channels and perturbative orders to the 
relative EW corrections and their sum. Thus, the blue line in plot (b) and the lower panel 
of plot (a) refer to the same quantity. Plot (c) shows the LO contributions of LL, RR, 
LR+RL processes and their sum; in the lower panel we display the relative EW corrections 
for each one of these classes. Panels (d), (e) and (f) respectively show the same kind of 
plot depicted in panels (b), but for the individual cases of LL, RR and LR+RL squark- 
antisquark production. 

Degenerate squark masses 

Slope Si: Fig. |] corresponds to the case of degenerate squarks, showing the dependence 
of squark-antisquark production on the common mass M g of all light flavor squarks. The 
inclusive cross section varies over two orders of magnitude in the considered mass range. 
The QCD corrections vary from 50% to 70% as M g varies from 1000 to 2500 GeV. The 
relative size of the total EW corrections increases with M g . and Fig. |3](b) clearly shows that 
the total amount of the EW corrections is the result of substantial cancellations among the 
different channels. This cancellation may be jeopardized by the photon PDF uncertainty, 
since the 57 -channel (with its large uncertainties) has a large impact on the total size of 
the EW corrections, as can be seen from the comparison of red and blue lines in the lower 
panel of Fig. ||](a). Looking at the different chirality combinations separately, in the case 
of LL production, mutual cancellations among the various channels keep the total EW 
corrections small in the entire M g range. On the contrary, in RR production the total 
EW corrections are positive and increase as M g increases. However, without the (possibly) 
large positive 57 -channel the EW corrections are negative in the entire considered mass 
range, and they reach — 20%(—5%) for M g = 2500 GeV for LL(RR) production. The total 
EW contributions to LR+RL production, which gives the largest part of the “inch” cross 

3 The physical masses are given by rriq L/R ~ l/r ^1 + O ^ M ™ z —^ . Thus, when e.g. A+^l/r = 
1/1.5/2 TeV, the physical mass and the soft mass parameters are equal at the 0.8/0.4/0.2% level. 
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section (see Fig. f|(c)), depend mildly on M ( j: they are negative and of the order of —5% in 
the entire range considered. As the ( 77 -channel does not contribute to LR+RL production 
these predictions are not affected by large photon PDF uncertainties. It is worth to notice 
that in LR+RL production the entire 0(a s a) contribution is in general zero, due to the 
different chiralities in the final state. Moreover, the size of the 0(a 2 ) contribution is 
negligible. Thus, in LR+RL production, the total EW corrections can be identified with 
the 0(a 2 a ) NLO EW corrections. 

Scan in M 2 : Results for the scan in M 2 are shown in Fig. [|. The QCD corrections 
to the inclusive cross section are independent of M 2 and of the order of 65%, while the 
total EW corrections are negative and small, of the order of 3%. The NLO EW corrections, 
0(a 2 a ), are negative and constant; the mild dependence on M 2 is induced by the EW tree- 
level induced contributions. As can be inferred from Figs. [|(c)-[|(f), the EW corrections 
do depend on the chirality of the produced squarks. In the case of LL production, they 
strongly depend on M 2 ranging from -3% to -10% as M 2 varies from 50 to 1200 GeV. This 
dependence is mainly induced by the 0 (a 2 ) contribution, which is suppressed for large 
values of M 2 . Both in LL and RR production the relative contribution of the 57 -channel is 
enhanced to 10%. At variance with LL production, the EW corrections in RR production 
are independent of M 2 , positive and of the order of 5%. As can be noticed in Fig. [|(f), for 
LR+RL production the only non-vanishing contribution to the total EW corrections is the 
0(a 2 a ) NLO EW, which is of the order of -5%. 

Scan in mg: The dependence of the QCD corrections on mg is mild and of the 
order of 60% in the entire range considered, as shown in Fig. ||. The EW corrections are 
small, negative and of the order of —4% in the low-m^ region. They increase to —2.5% for 
nig — 2000 GeV. As can be inferred from Figs. ||(b), the increment is mainly due to the 
positive yield from the 0(a 2 ) contribution and the 57 -channel. Again, the 57 -channel and 
the corresponding PDF uncertainty are of the same order as the EW corrections themselves. 
As shown in Figs. [|(c) and [|(d), the EW corrections are more important in the case of LL 
production. In the large mg region, the 0(a 2 ) and the 57 -channel contributions become 
dominant and render the EW corrections positive and small, of the order of 1%. As 
shown in Fig. 1(e) and 1 (f) the total EW corrections to RR and LR+RL production are 
respectively of the order of —2% and —5% in the entire region considered. Again, in LL and 
RR production the photon PDF could substantially alter the size of the EW corrections, 
while for LR+RL production the photon PDF does not contribute. 

Non-degenerate squark masses 

In the following scenarios the left- and right-handed soft squark masses rriq R and rriq L are 
treated as non-degenerate parameters. Here, NLO QCD corrections cannot be computed 
by Prospino and are not included in the analysis. 

Slope S 2 : For this scan, displayed in Fig. 1 we set rriq R = 2niq L ; hence, RR production 
is negligible and the cross section for LR+RL is much smaller than for LL production. 
Therefore EW corrections to the inclusive squark-antisquark production are qualitatively 
very similar to the EW corrections to LL production in the slope S\. The overall size of 


- 14 


the EW corrections is the result of mutual cancellations among the different channels and 
their dependence on Mq is mostly determined by the 57 -channel. Thus, the photon PDF 
uncertainty may substantially alter these corrections. RR production is strongly suppressed 
by the high mass of the produced squarks, i.e. the scan covers values of rriq R in the range 
[2 TeV, 5 TeV]. In this regime, the photon PDF uncertainty renders the EW corrections 
to RR production unreliable. The clear sign of their unphysical behavior is the “bump” at 
Mq = rriq R /2 = 1.7 TeV, which is due to the photon PDF and not from matrix elements. 

Slope S 3 : In this scenarios (see Fig. j|) with rriq L = 2 rriq R , both LO predictions 
and higher-order corrections are mostly determined from RR production. The summed 
EW corrections are positive and increase as Mq increases, reaching 10% in correspondence 
to Mq — 2000 GeV. The mass hierarchy is inverted wrt. the one of the slope S 2 , so 
the qualitative discussion for the LL(RR) contribution in S 2 applies here for the RR(LL) 
contribution. The main difference with S 2 , as can be seen from Fig. |7|(e), is the fact that 
there are no large cancellations among the different EW contributions for RR(inclusive) 
production. 

Slope S 4 : In this scan (Fig. ||) the mass of the right-handed squark mq R is kept 
fix at rriq R = 1.5 TeV, while rriq L = Mq is varied. As can be seen in the lower panel of 
Fig. |H(c), the hierarchy of the LO predictions for LL, RR and LR+RL depends on Mq. 
Consequentially, also the EW corrections and the individual perturbative orders receive 
the dominant contributions from LL, RR or LR+RL depending on the value of Mq. In RR 
production, Fig. [ 8 |(e), the EW corrections are constant by construction and corresponds to 
those of the spectrum in Table (2| Also the summed LR+RL production, Fig. |(f) , does 
not show a visible dependence on Mq. LL production shows a very similar behavior as 
observed for the slopes Si and £ 2 - 

Slope S 5 : In this scan, referring to Fig. ||, the mass of the right-handed squark rriq L 
is kept fix at rriq L = 1.5 TeV, while rriq R = Mq is varied. Thus the LO cross sections are 
those of Slope £4 with LL and RR exchanged, cfr. Fig. ||(c) and Fig. ||(c). The qualita¬ 
tive behavior of the electroweak corrections to the different production channels can be 
understood as an interchange of LL with RR with respect to Slope £4. 

3.3 Differential distributions 

In the previous subsection we studied the numerical impact of the EW corrections to the 
total cross section for squark-antisquark production at the LHC. It is well known that, 
at high energies, Sudakov-type logarithms can enhance the EW contributions. Therefore, 
in the following we study the impact of the EW corrections on three kinematic distribu¬ 
tions: the transverse momentum of the produced squark, the maximal pseudo-rapidity of 
the squark and the anti-squark, and the invariant mass of the squark-antisquark system, 
defined as 

\ Vg a if 

Pt = Pt q a ; r]=\ 

y if 

respectively. The quantities pj, pxj and i]j denote the four-momentum, the transverse 
momentum and the pseudo-rapidity of the particle j. It is worth to remind that these 


k«l > \V? p \ 

kj < k; I 


Mu 


- \J (pq a +Pq%) , ( 3 - 4 ) 
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distributions cannot be directly observed at the experimental level; a realistic phenomeno¬ 
logical evaluation requires the combination of the production process with the decays of the 
squarks. Nevertheless, it is important and useful to identify the kinematic dependence of 
EW corrections at the production level, before performing a complete simulation including 
squark decays and acceptance cuts. 

For illustration of the effects, we consider, as a benchmark, the scenario defined in 
Table ||, but now with M^l = = 2000 GeV, i.e. a scenario within the reach of 

Run-II of the LHC. The numerical results for distributions in px, i] and M\ m are collected 
in Fig. 10, Fig. [II] and Fig. |l^ respectively. In the left plots we always show the absolute 
predictions at LO with and without the EW corrections. The right plots contain the 
relative size of the EW corrections with respect to the LO prediction, and the breakdown 
to their various individual contributions. From top to bottom, the plots are arranged to 
display the predictions for inclusive, LL, RR, and LR+RL production. The overall behavior 
with respect to the different kinematic observables is strongly dependent on the chirality 
configuration of the produced squarks. 

The EW corrections to the px distribution, shown in Fig. |l(], are largest for the LL 
squark configurations, owing to their hypercharge. In particular they are negative and reach 
20% and more for large px- The LR+RL configurations, which dominate the inclusive cross 
section, are influenced by —10% for large pt■ In the RR configurations, the EW corrections 
always stay below 5%. 

The 0(a s a ) and 0(a 2 ) Born contributions, relevant only for LL and RR production, 
have always opposite sign and similar magnitude, resulting in mutual cancellations in most 
of the p-T-spectrum. The photon-induced production channel, again only relevant for LL 
and RR production, can yield sizable contributions both at small and at very large pt■ In 
particular, in the case of RR production it tends to overcompensate the negative NLO EW 
contribution. 

Turning now to the EW corrections in the pseudo-rapidity p of the produced squarks 
(Fig. the origin of the large photon-induced contributions at small px becomes evident. 
The photon induced production modes show a different angular dependence as compared 
to the LO prediction or the NLO EW, i.e. they are strongly enhanced in the forward region 
at large |? 7 |, while the NLO EW corrections are mildly enhanced for small \p\. Due to this 
large dependence on r/, the relative corrections induced by photon-induced contributions 
on experimental rates can strongly be affected by the details of the analysis cuts on the 
squark decay products. 

Finally we turn to the distributions in the invariant mass of the produced squark- 
antisquark pair (Fig. 12). Again we observe a partial mutual cancellation between the 
different contributions to the EW corrections. Moreover, the NLO EW corrections are not 
significantly enhanced at large invariant masses. 


4. Conclusions 

In this paper we have presented the first phenomenological study of the NLO EW correc¬ 
tions to squark-antisquark production at the LHC including all production channels and 
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chirality combinations of the produced squarks. 

Our analysis has shown that electroweak contributions to q a q'p production are not 
negligible even at the inclusive level. The 0(a 2 a) NLO EW contributions are negative, in 
general sizable and increase in relative size with the mass of the produced squarks. For 
the production of a left-handed squark and antisquark pair they reach —20% for squark 
masses of 2 TeV. For the production of two left- or two right-handed squarks these negative 
contributions are typically (over-)compensated by large tree-level contributions from the 
photon-gluon initial state, which also increase in relative size for large masses. However, 
these contributions are accompanied by very large intrinsic PDF uncertainties, which may 
substantially alter the size of the electroweak corrections and, for large masses, the accuracy 
of the total predictions. The dependence on the remaining MSSM parameters is found 
to be very weak for the NLO EW 0(a 2 a) and the photon-induced contributions. The 
subleading Born contributions of 0(a s a ) and 0(a 2 ) are not negligible and show a moderate 
dependence on mg and M 2 . 

Besides inclusive cross sections, the electroweak corrections to squark-antisquark pro¬ 
duction are also investigated at the differential level. The 0(a 2 a) NLO EW corrections 
show a typical Sudakov behavior, i.e. increasing influence at large transverse momenta of 
the produced squarks. As on the inclusive level, these corrections are partly compensated 
by large contributions from photon-induced production. The origin of these large contribu¬ 
tions can clearly be attributed to a qualitatively different behavior of the photon-induced 
channels at large pseudo-rapidities. 

In general, in studies for squark production as well as in corresponding data analyses, 
when EW effects are usually neglected for reasons of simplicity, the size of the EW contri¬ 
butions calculated and visualized in this paper, can serve as an estimate of the uncertainty 
of the theoretical predictions, on top of the uncertainty resulting from the QCD side. 
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Figure 3: Scan over Mg, which corresponds to Slope Si. The value of the other parameters 
are collected in Table |2| 
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Figure 4: Scan over M 2 . The value of the other parameters are collected in Table 
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Figure 6: Slope S 2 ' scan over Mq. The value of the parameters not involved in the scan 
are collected in Table |2| 
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Figure 7: Slope S 3 : scan over Mq. The value of the parameters not involved in the scan 
are collected in Table |2| 
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Figure 8: Slope S 4 : scan over Mq. The value of the parameters not involved in the scan 
are collected in Table |2| 
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Figure 9: Slope S 5 : scan over Mq. The value of the parameters not involved in the scan 
are collected in Table |2| 
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Figure 10: Differential distributions in the transverse momentum pj- of the produced 
squark. 
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Figure 11: Differential distributions in the maximal pseudorapidity q of the produced 
squarks, as defined in eq. (3.4). 
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Figure 12: Differential distributions in the invariant mass M- im of the produced squark- 
anti-squark pair. 
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